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ABSTRACT 



In the advent of next generation gamma-ray missions, we present general prop- 
erties of spectral features of high-energy emission above 1 MeV expected for a class 
of millisecond, low magnetic field (~ 10 9 G) pulsars. We extend polar-cap model cal- 
culations of Rudak & Dyks (1999) by including inverse Compton scattering events in 
ambient field of thermal X-ray photons and by allowing for two models of particle 
acceleration. In the range between 1 MeV and a few hundred GeV the main spectral 
component is due to curvature radiation of primary particles. Synchrotron compo- 
nent due to secondary pairs becomes dominant only below 1 MeV. The slope of the 
curvature radiation spectrum in the energy range from 100 MeV to 10 GeV strongly 
depends on the model of longitudinal acceleration, whereas below ~ 100 MeV all 
slopes converge to a unique value of 4/3 (in a vT v convention). The thermal soft X- 
ray photons, which come either from the polar cap or from the surface, are Compton 
upscattered to a domain of the VHE and form a separate spectral component peaking 
at ~ 1 TeV. We discuss observability of millisecond pulsars by future high energy 
instruments and present two rankings relevant for GLAST and MAGIC. We point to 
the pulsar J0437-4715 as a promising candidate for observations. 

Key words: gamma-rays: theory, observations - pulsars: general 



1 INTRODUCTION 

Millisecond pulsars are thought to be uninteresting as tar- 
gets for gamma-ray experiments. This opinion partly draws 
from an argument that typical values of B (or rather their 
dipolar components) inferred from P and P are in the range 
of 10 s to 10 9 G, in contrast with the range of 10 11 to 10 13 G 
for classical pulsars. If gamma-ray radiation is to be a man- 
ifestation of magnetospheric activity, it seems quite natu- 
ral to think of millisecond pulsars in terms of B as of a 
scaled-down version of classical pulsars, which by themselves 
are weak gamma-ray s ources; only seven of th em have been 
firmly detected so far (Thompson et al., 1997). Unsuccesful 



there are 9 millisecond objects (Becker & Triimper, 199!:) 



observational campaigns in the past, mainl y by CGRO in- 
struments e.g. EGRET - (Fierro et al., 1995), confirmed this 
line of thinking. J0437-4715, with one of the highest spin- 
down fluxes among all pulsars, thus a very promising ta rget 



wasn't even included in the COMPTEL priority list (Car- 
raminana et al. 



1991 



However, the 'scaling-down' argument becomes unjusti- 
fied by the results of X-ray observations. Within a sample of 
some 30 pulsars detected with ROSAT and other satellites 



with apparent X-ray properties similar to classical pulsars. 
In particular, the millisecond pulsars ob ey th e same aston- 
ishing relation of Becker & Triimper (1997) between the 
spin-down luminosity L sc \ and the soft X-ray luminosity 
L x that classical pulsars do. In most cases the X-ray spectra 
are presumed to be non-thermal in character, but whether 
they represent low-energy tails of putative gamma-ray emis- 
sion is not known. 

Rather few papers have been dedicated to the question 
of gamma-rays in the c ontext of millisecond pulsars on theo- 
retical side. Wei et al. ( 1996 ) presented broad-band spectra 
of both pulsed and unpulsed emission expecte d in t heir ver- 
sion of outer-gap model. Sturner & Dermer (L994) offered 
predictions for gamma-ray luminosity L 7 in some popular 
models, but that had been done just by scaling-down the 
values of B in simple analytical formulae for L 7 derived for 
classical pulsars . In t he framework of polar-cap scenarios, 
Rudak & Dyks (1998) pro posed a modification of the model 
by Daugherty & Harding (J198§). This modification incorpo- 



rates a contribution of electron-positron pairs to L T , which 
becomes a non-monotonic function of B (see Figures 7 and 
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8 of ( Dyks, 1998 )). One of the consequences of this modifi- 
cation is that expected values of L 1 for millisecond pulsars, 
including J0437-4715, are typically about one or two orders 
of magnitude lower than for classical pulsars. Quantitatively 
similar estimate, though for differ ent reason, comes from a 
model of Dermer & Sturner ( 1994 ). 

However, if the estimates above are correct then some 
millisecond pulsars, including J0437-4715, should be de- 
tected with next-generation experiments, providing thus a 
testing ground for magnetospheric models relevant for mil- 
lisecond pulsars. The successor of EGRET - GLAST - 
with a possible launch in the year 2005, is expected to 
have sensitivity above 100 MeV about 30 times higher then 
EGRET , and it will reach high-energy limit at 300 GeV 
- closing thus for the first time the energy gap between 
the VHE (very high energy) domain accessible with ground- 
based Cherenkov techniques a nd the HE (high ene rgy) do- 
main of satellite experiments (Kamae et al., 1999). More- 
over, since we expect the millisecond pulsars to radiate the 
electromagnetic energy mostly in the range between 10 GeV 
and 100 GeV, this class of objects sho uld be of interest fo r 
the upcoming IAC telescope MAGIC (Blanch et al., 1998). 

Our aim is to present major spectral features in the en- 
ergy range above ~ 10 MeV of pulsed gamma-ray emission 
which arises due to both ultrarelativistic primary particles 
(electrons) and cascades of secondary particles above po- 
lar caps of millisecond pulsars with dipolar magnetic fields 
of ~ 10 9 G. In particular, we allow for different structures 
of the accelerating potential, and we show how this affects 
the slope of the predicted energy spectrum in the range of 
100 MeV to 100 GeV. Also, Compton upscattering of am- 
bient X-ray photons off ultrarelativistic electrons is taken 
into account. In Section 2 we present details of modelling 
the processes responsible for the formation of gamma-rays 
in millisecond pulsars. Section 3 describes spectral proper- 
ties of two distinct components, due to curvature radiation 
(CR) and Compton scattering (ICS). Conclusions with a 
discussion of J0437-4715 come in Section 4. 



2 THE MODEL OF HIGH-ENERGY 
RADIATION 

We use a polar cap model, with most ingredients as postu- 
lated by Daugherty & Harding (1982) in context of classical 
pulsars. High-energy emission is a superposition of curvature 
radiation (CR) of ultrarelativistic primary electrons, and 
synchrotron radiation (SR) of secondary particles (e* pairs) 
created via magnetic absorption of photons. Additional fea- 
tures include: 1) two models of electron acceleration, 2) in- 
verse Compton scattering (ICS) of thermal X-ray photons 
on electrons. The first feature is important for the history of 
cooling of primary electrons and consequently it affects the 
slope of gamma-ray spectrum above 100 MeV. The second 
one (ICS) is not important for overall cooling of electrons, 
but the resulting TeV component in the gamma-ray spectra 
is of potential interest. A field of ambient photons necessary 
for the ICS is likely to be present within the magnetosphere 
of some millisecond pulsars, according to some interpreta- 
tions of available X-ray and EUV observations (e.g. Zavlin 
& Pavlov 1998, Edelstein et al. 1995). We list the details of 
the model in the following subsections. 



2.1 Geometry 

We assume that the geometry of the magnetic field of a neu- 
tron star is well described by a static, axisymmetric dipole 
In polar coordinates its field lines satisfy 

-2 



rsin 6 — Ro, 



(1) 



where the dipole constant 7?o defines a set of field lines 
which differ by azimuthal angle only. Polar caps are defined 
as (two) regions on the neutron star surface crossed by all 
lines for which the condition Ro > Ri c is satisfied, where 
R\ c = cP/2n is a light-cylinder radius for a spin period P. 
The polar cap radius is then r pc = Rns x (i? n s/fiic) 1//2 - where 
R na is the radius of the neutron star for which we assume 
the canonical value of 10 6 cm. Beam particles (primary elec- 
trons) are ejected from the outer rim of the polar cap and 
move along the magnetic field lines. 

When modelling an extended source of soft photons, 
which are to be targets for electrons in Compton scatter- 
ings, we consider two cases of its geometry. In the first case, 
the soft photons come from the hot polar cap with a uni- 
form temperature T pc ii 10 6 K. In the second case the soft 
photons come from the entire surface of the neutron star, 
however now the temperature is lower: T sur f JSJ 10 s K. 



2.2 Electron acceleration 

Primary electrons are injected along magnetic field lines into 
the magnetosphere from the outer rim of the polar cap. The 
structure of longitudinal electric field responsible for their 
acceleration remains unknown. However, the electric field 
must be strong enough to ensure electron positron pair cre- 
ation. The presence of pairs is a necessery condition for a 
magnetized neutron star to be a radio pulsar. Thus the re- 
quirement of creating pairs in the magnetosphere is the pri- 
mary condition that the model has to satisfy. 

We use two simple models for particle acceleration. The 
first one (denoted model A) assumes instant acceleration to 
ultrarelativistic energy _EJ, nlt at the moment of electron in- 
jection on the surface of the neutron star. The values of E e nlt 
are set to be just above the threshold value for creation of 
secondary pairs _B™ m which in turn is a function of mag- 
netic field strength at the polar cap B pc and the spin period 
P (for details see Rudak & Dyks 1998, and Dyks 1998). 
For a pulsar with B pc — 10 9 G and P = 3ms, the initial 
energy is E [ e nit = 1.07 x 10 7 MeV. For J0437-4715 we took 
B pc = 7.4 x 10 8 G, P = 5.75ms, and £^ nit = 1.54 x 10 7 MeV, 
where we used the magnetic field derived from P kinemati- 
cally corrected by Camilo et al. (1994). 

For comparison we consider another model (model B) of 
particle acceleration, in which the longitudinal electric field 
£ is set at a constant level with an exponential cutoff: 



- V h 
b = exp 



(2) 



where the characteristic scale-height of the electric field is 
assumed to be equal to the radius of the polar cap r pc . For 
a pulsar with B pc = 10 9 G and P = 3ms the value of £ 
at the polar cap surface was set at the level of Vo/r pc ~ 
10 9 V cm -1 to ensure a total yield of secondary pairs per 
primary electron similar to that in model (A). The initial 
value of the electron energy E e u,t = 10 MeV was assumed. 
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Figure 1. Cumulative optical depth r e for non-magnetic scattering of an electron injected at the surface level (h = 0) and crossing a 
blackbody radiation field. The left panel is for the radiation fields due to hot polar cap (T cap = 10 6 K and 3 X 10 6 K). Solid lines are 
for instant acceleration with electron energy E\ n ^ t = 10 7 MeV. Dashed lines are for continuous acceleration in electric field £ (eq.^). The 
right panel is for the radiation fields due to hot stellar surface (T sur f = 2 X 10 5 K and 5 X 10 5 K). The pulsar period is P = 3 X 10 _3 s 
and we used 10 5 electrons in the simulation. 



As electrons accelerate they move along the magnetic 
field lines and are slowed down by the curvature radiation 
reaction. We calculate the electron energy budget as a func- 
tion of the distance travelled in the magnetopshere. Only a 
small fraction (less than one percent, see Figure |l]) of the 
electrons will scatter on thermal photons originating on the 
polar cap or on the stellar surface. Thus scattering is not an 
important electron deceleration process. 



2.3 Curvature and synchroton processes 



We follow the model described by Rudak & Dyks (199£ ) 



who calculated broad band energy spectra of high-energy 
emission for two distinct pulsar groups: low-magnetic field 
pulsars ('millisecond pulsars') and classical pulsars. 

The dominant cooling process for primary electrons is 
the curvature radiation: 



(7e)cr = (2/3)7*P 2 CT , 



(3) 



where 7 e is the electron Lorentz factor, p is the curvature 
radius, and ro is the electron radius. When considering sec- 
ondary electrons (pairs) one has to consider also the syn- 
chrotron radiation. The rate of SR cooling is 



/ • \ 2 r\ 2 • 2 , 2 

(7e) a r = ~- B Sill V7e , 

6 m e c 



(4) 



where ip is the pitch angle of the pairs. The photon spectra 
due to curvature and synchrotr on pr ocesses are calculated 
following Daugherty & Harding (1982). A detailed dicussion 
of the role of these processes in the formation of ga mma-ray 
spectra of pulsars is given in Rudak & Dyks (1999). 




10 4 10 5 
E[MeV] 

Figure 2. The dotted line shows the intrinsic spectrum of the 
inverse Compton photons i.e. with transfer effects neglected. The 
solid line shows the spectrum of magnetically reprocessed pho- 
tons as seen by an external observer. Here we assumed: thermal 



emission from a polar cap with the temperature T ci 



: 3xl0 6 K, 



surface magnetic field B = 10 9 G, pulsar period P = 3 X 10 
and instantaneous electron acceleration. 



2.4 Scattering 

Our numerical treatment of scattering process follows 
Daugherty & Harding (1989). As the electrons travel up in 
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T sur[ = 5x10 5 K 



10 3 10 3 10 4 10 5 10" 10 7 10" 

h[cm] 




10- 6 r 
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Figure 3. The thick solid lines in the left panels show the dependence of the mean energy of unscattered primary electrons as a function 
of the height. The top left panel corresponds to the case of instantaneous acceleration (model A), while the bottom left panel corresponds 
to case of accleration described by eq. (j^J) (model B). The dots show the energies of the Compton upscattered photons and height where 
the scatterings took place (in these simulations we propagated 10 5 primary electrons in the magnetosphere). The thin solid line shows 
the optical depth to magnetic absorption equal to unity, according to equation (|ll|). The right panels show the energy spectra formed 
for model A (top right panel) and model B (bottom right panel). The soft photons originate on the surface of the neutron star with 
the temperature T sur f = 5 X 10 5 K. The CR and the ICS spectral components are due to a single primary electron injected at the 
outer rim of the polar cap of the standard millisecond pulsar. We normalise the spectra with iJ™ ax - maximum energy of the electron, 
mediated between acceleration and cooling rates. For the top left panel, £™ ax is equal to the initial energy Ei nit = 1.07 x 10 7 MeV. 
Pulsar parameters are B pc = 10 9 G, and P = 3ms and £^ nit = 1.07 X 10 7 MeV for model A. 



the magnetosphere they may scatter off thermal photons 
originating from the stellar surface or the polar cap. For 
millisecond pulsars the magnetic field at the surface is typ- 
ically 10 9 G or less, i.e. much weaker then the critical field 
B cr it — 4.414 x 10 13 G. Thus the scattering process will be 
well described by the Klein-Nishina relativistic nonmagnetic 
cross section: 



&KN — 



3(TT 

4a; 



■ Ax - 8 



ln(l + x)- 



1 8 _ 1 
+ 2 + x 2(1 +x) 2 



(5) 



where: x = 2p v k v , p" is the electron four-momentum, and 
k v is the photon four-wave vector. The scattering rate (in 
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the "lab frame" K) by an electron travelling with velocity 
(3 in a photon field is given by 



1Z = c / dQ I £ 2 ciecrA']v/e(l — P ■ hi) 



(6) 



where Q is the solid angle subtended by the source of soft 
photons, (tkn is the total cross section, f e is the photon 
occupation density, and hi is the unit vector (hi = k/\k\) 
representing the direction of the incoming soft photon. Tak- 
ing the Planck function for e 2 / e this equation reads: 



n = 



c /m e c\ : 
4^3 \~) 



dQ I di okn{^ — fi-hi) 



exp e/T - 1 



(7) 



where e is the photon energy in the units of electron mass, 
T is the photon temperature in the units of electron mass. 

We draw the energy and direction of the incoming pho- 
ton from the distribution of e 2 / e , and Q, respectively (as 
used in eq.6). Then we transfer them to the electron in- 
stantaneous rest frame K e and there we find the photon 
scattering angle using the Klein-Nishina differential cross 
section 



do j. 



3ot 



/' 



dfj,d(f> 



167T [ 1 + f(1 _ M) ] 2 



4(1 + /i*) [1+ |(l - M )] 



(8) 



where fi is the cosine of 9. The energy e of the scattered 
photon in K e is obtained from the Compton formula 

(9) 



1 + 



r(l-M) 



Finally, the energy of the outgoing photon is transformed 
back to K. 

At each step, as the location of the electron changes, 
the solid angle Q is corrected for the finite size of the source 
of soft photons (the polar cap or the stellar surface). Thus, 
we take into account the effect of the decreasing photon 
number density with height (geometrical dilution) and then 
we calculate the scattering rate 1Z at a given height above 
the neutron star surface. 



2.5 Magnetic photon absorption 

Despite low magnetic field values, two conditions for the 
process of pair creation via magnetic photon absorption (7 + 



B 



the energy threshold condition, and the high 



optical thickness r 7 s of the magnetosphere - may well be 
met in the context of millisecond pulsars. For instance, in 
the model spectra calculated numerically by Rudak & Dyks 
(1999) photons with the energy ~ 10 2 GeV or higher were 
subject to the magnetic-absorption reprocessing. 

In our simulations we use the absorption coefficient as 
described by Erber (1966): 



1 a Bj 



2 A c B c 



-T(X) 



(10) 



where a is a fine structure constant, A c is a Compton wave- 
length, B± is the component of the magnetic field perpen- 
dicular to the photon momentum, and x 



2 B crit m e c' 2 

the Erber parameter \- F° r the function T(x) we use its 



approximation T(x) ~ 0.46 exp (— 4//3x), valid for x ^ 0.2 
(for x 0.2 this approximation starts to overestimate rj). 
Th e function f is the near-thre shold correction introduced 
by ( Daugherty fe Harding, 198S ) . In the case of millisecond 



pulsars equation (|10|) is a good approximation even with 
/ = 1, since magnetic photon absorption occurs well above 
the threshold. Electron-positron pairs created through the 
magnetic absorption radiate then via the synchrotron pro- 
cess. 

For demonstrative purposes it is useful to present the 
condition r 7 _g > 1 in an analytical way. It is straightforward 
to show that a photon created with a momentum parallel to 
the local magnetic field line at a height h above the neutron 
star surface will undergo magnetic absorption if its energy 
satisfies approximately the following inequality: 



E > 10' 



uo-W 

\ Rns ^ 



Bn 



U0 9 G 



10 6 cm 



-1/2 



5/2 



MeV. 



(11) 



This formula is valid for dipolar magnetic field, and the field 
line is attached to the outer rim of the polar cap, and it 
is in very good agreement with numerical calculations of 
magnetic photon absorption. 



3 RESULTS 

Electrons are either accelerated instantly (model A) or ac- 
celerated by the electric field given by equation | (model 
B), and decelerated mainly by the curvature radiation. In 
model B the deceleration rate is insignificant near the polar 
cap because it increases with the electron energy and E 1 ™ 
is small. On the other hand, the acceleration rate does not 
depend on the electron energy and it ceases at the height 
h ~ r pc (eq.^) . There exists a region in between where both 
competing processes balance. An example calculation of the 
electron energy is shown in Figures || and ^. The top left 
panel of each figure presents electron energy within model 
A, with Ei nit = 1.07 x 10 r MeV. The bottom left panel shows 
the electron energy within model B; the electron reaches a 
maximum energy of a few TeV at the height h ~ 4 x 10 4 cm 
(i.e. h/r pc ~ 1/6). 

The optical depth r e for nonmagnetic scattering of an 
electron is small (Figure |l|) and thus scattering is not a rele- 
vant decelaration process of electrons. The optical depth for 
scattering scales as oc T 2 , where T is the temperature of the 
radiation field. There are several factors that influence the 
scattering rate as a function of the height above the surface 
of the neutron star h. The photon density decreases with 
height, but the electron spends much more time in the up- 
per part of the magnetosphere, thus increasing the chance of 
scattering there. We find that the latter effect is dominant 
and most of the scatterings take place in the upper part 
of the magnetosphere between 10 s and 10 7 cm, where the 
electron energy is close to the value of a few TeV. This is il- 
lustrated in Figures ^| and ^ (left panels) , where the electron 
energy is plotted as a function of h. We show the energies of 
the ICS photons and the locations of the scattering events. 
The intrinsic spectrum of the ICS photons has a width of 
about one decade in energy centered around a few TeV, cor- 
responding to the range of electron energies in the region 
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Figure 4. Same as Figure pi but for the case of thermal emission from the polar cap with the temperature 3 X 10 6 K. 



where most of the scattering occur. In Figure g we illustrate 
the effects of reprocessing in the pulsar magnetosphere (tak- 
ing model A as an example) and present both, the intrinsic 
spectrum of ICS photons (dotted line) and the spectrum 
after reprocessing due to magnetic absorption (solid line). 
Magnetic absorption removes the high energy part of the 
intrinsic ICS spectral peak, decreasing the strength of the 
peak by a factor of more than 10, and redistributing the 
ICS photon energy to lower part of the spectrum. The ICS 
photons subject to magnetic absorption reprocessing are lo- 
cated in Figures ^| and ^ (left panels) above the thin line 
corresponding to r 7 _g = 1 (equation ||ll|). In the case when 
the soft photon source is the hot polar cap most of the ICS 
photons will be reprocessed; see Figure kl However, in the 
case when the photon source is a warm surface of a neutron 
star most of the scattering events take place at a distance 
of a few neutron star radii from the surface and most of the 



ICS photons will escape freely; see Figure |§| The escaping 
ICS photons have typical energies equal to the electron en- 
ergy at the height of few neutron star radii above the stellar 
surface and thus the spectrum of the ICS photons exhibits 
a strong peak at the energy of ~ 1 TeV. The height of the 
ICS peak depends on the number of the soft thermal pho- 
tons in the magnetosphere, and so it is directly related to 
the temperature and the radiating area (polar cap or entire 
stellar surface). 

Another factor that influences the scattering rate is the 
soft photon flux anisotropy. Such anisotropy results from 
opacity angular dependence in strong magnetic field, and 
cannot be neglected in the case of classical pulsars with 
B ~ 10 12 Gauss. Detailed models of magnetized atmo- 
spheres show that in the case 10 9 Gauss fields the so ft pho- 
ton flux is almost constant up to about 60 degrees (Pavlov 
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Figure 5. Millisecond pulsars ranked according to their gamma 
ray fluxes above 100 MeV expected within model A. Each verti- 
cal bar spans the range in the flux when the energy of injected 
primary electrons i?g mt varies by a factor of two. The horizontal 
dashed line represents the GLAST sensitivity. 



et al., 1994). This result justifies our simplified treatment of 
the soft photon field as being isotropic. 

The dominant spectral component in the energy range 
from 10 MeV out to its cutoff at <; 100 GeV is due to 
curvature radiation (see the energy spectra in Figures ^ 
and ^j, right panels) . This cutoff energy is significantly higher 
than that expected in some classical pulsars as well as me a- 
sured in some of them (ss 10 GeV, Thompson et al. ( 1997 )). 
This is because of smaller curvature radii implying higher 
values of the characteristic photon energy, and also much 



weaker magnetic absorption {■yB — ► e ). The slope a of ra- 
diation energy spectrum per logarithmic energy bandwidth 
(E 2 dN/dE oc E a ) below e brcak « 100 MeV is insensitive 
to model of acceleration, and is close to 4/3 (see Rudak & 
Dyks 1999 for discussion of eb rea k)- Above ~ 100 MeV the 
spectrum assumes a slope which depends on the acceleration 
model. For model A the slope between 100 MeV and 10 GeV 
is 1/3, in agreement with simple analytical estimates for mo- 
noenergetic source function of beam particles. The slope for 
model B was found numerically to be equal 0.84 (Bulik & 
Rudak, | 1999| ). The electron acceleration in model B makes 
the electron spectrum harder in comparison to the case of 
model A, which in consequence leads to a harder CR spec- 
trum. Note, that the total energy contained in gamma rays 
is larger in model B than in model A (compare right-hand 
panels of either Figure |^ or Figure . In model A the initial 
energy of electrons is large enough to enable pair creation, 
while in model B they need to overcome the CR losses to 
attain £" lax . The work done by the electric field on the elec- 
tron in model B is about 20 times larger than _E™ ax . 



Figure 6. Similar to Fig. 6 but the ranking is for gamma ray 
fluxes above 100 GeV, and the horizontal dashed line represents 
the MAGIC sensitivity. 



4 DISCUSSION 

We calculated energy spectra of pulsed HE and VHE emis- 
sion from millisecond pulsars for two models of beam par- 
ticles acceleration. The spectra are a superposition of cur- 
vature (CR), synchrotron (SR), and Compton upscattering 
(ICS) components. The curvature component dominates the 
region between 1 MeV and 100 GeV. The slope of the spec- 
trum in the range between 100 MeV and 10 GeV is sen- 
sitive to the details of the electron acceleration process. 
The synchrotron component becomes important only below 
~ 1 MeV (Rudak & Dyks 1999). The ICS component has a 
narrow peak around 1 TeV and its strength is related to the 
number density of the soft photons in the magnetosphere. 
Two sources of soft photons were considered: the hot polar 
cap with a uniform temperature T pc ?t 10 6 K, and the en- 
tire surface of the neutron star, with a lower temperature: 
T surf £ 10 5 K. 

Several new g amma-ray tel escopes will be operational 
in the near future (Weekes, 1999) and thus it is interesting to 
discuss the observability of the HE radiation from millisec- 
ond pulsars. For our ranking we have selected the millisec- 
ond pulsars from the pulsar catalog of Taylor et al. (1995) 



using the following criteria P < 60 ms and P < 10 _1 \ For 
each pulsar the expected integral fluxes above 100 MeV and 
above 100 GeV was then calculated. Model A had been cho- 
sen for these calculations as a more conservative one. More- 
over, two cases were considered for each pulsar: in addition 
to the case with the energy of injected primary electrons 
E™ 1% being chosen according to section 2.2 (see also eq.8 of 
Rudak & Dyks 1998), we also allowed the electron energy 
to be two times higher. This was done to present how sen- 
sitive the results are to our choice of th e elec tron's energy. 
We used the results of Rudak & Dyks (1998) to callibrate 
the fluxes of primary electrons (note that the spectra pre- 
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Figure 7. Integral spectral fluxes of photons calculated for 
J0437-4715 at the distance of 140 pc. Model A is used with two 
following values of £j, nit : 1) 1.54 X 10 7 MeV (after Rudak & 
Dyks 1998, see also section 2.2), 2) a value two times higher, i.e. 
3.08 X 10 7 MeV. These two values give the spectra (solid lines) 
which limit the shaded region from the bottom and the top, re- 
spectively. The soft photons used for ICS come from the entire 



stell ar su rface with the temperature 4 X 10 5 K (Edelstein et al 



1995). Sensitivities of several high energy detectors are also indi- 
cated. 



sented in Figures ^ and ^, are per single electron) . Ten best 
candidates in each band are presented in Figures |E| and ^ 
respectively. The first two candidates for GLAST observa- 
tions (Figure ||), B1534+12 and B1821-24, are objects with 
relatively high magnetic fields, which place them in a transi- 
tion region between the domains of classical and millisecond 
pulsars. Because of the magnetic field strength and conse- 
quently high magnetic absorption in these two objects, the 
high energy cutoff lies below 100 GeV and therefore they 
disappear from the ranking of Figure ^. 

The third best candidate for GLAST observations is 
J0437-4715, a pulsar located very close to the Earth at the 
distance of 140 pc, with the spin period P — 5.75 ms, and 
the magnetic field at the polar cap B ^ c — 7.4 x 10 s G in - 
ferred from kinematically corrected P flCamilo et al., 1994|) . 
J0437-4715 is also our primary candidate for observations 
in the range above 100 GeV with the future Cherenkov 
telescopes. The integral photon spectra of J0437-4715 ob- 
tained in the course of calculations for the purposes of the 
ranking are presented in Figure [?] along with the sensitiv- 
ities of GLAST, CANGAROO, as well as MAGIC in its 
'Large Zenith Angle' mode. The expected photon flux in 
the EGRET energy range (i.e. between 0.1 and 10 GeV) is 
between w 2.0 x 10" 



,andRi 2xl0~ 8 cm~ 2 s~ 1 (for 
E [ e nit = 1.54 x 10 7 MeV and El nit = 3.08 x 10 7 MeV, respec- 
tively) and falls below the upper lim it of 2.1 x 10~ 7 cm " 2 s _1 
placed with EGRET for this object ( [Fierro et al., 1995| ). Fig- 
ure shows that the flux expected from J0437-4715 in the 



GLAST energy range (0.1 to 2 00 GeV) is above its sensitiv- 
ity limit of 3 x 10" 9 cm" 2 s" 1 flKamae et al., 199S| ). 

In the list of potential targets for atmospheric 
Cherenkov telescopes, presented in Figure |(| J0437-4715 
clearly stands out as the strongest source, mainly because 
of its proximity. The next source B 1257+ 12 is also rather 
strong and should be detected by the telescopes in the north- 
ern hemisphere like MAGIC and VERITAS, and also it could 
be seen by CANGAROO-III. The remaining objects have 
the expected fluxes between a few times 10 -13 and a few 
times 10 _12 cm _2 s _1 , which places them in a class of possi- 
ble but uncertain targets. 

The 100 GeV part of J0437-4715 spectrum can be 
probed by Cherenkov telescopes in the southern hemisphere 
like CANGAROO-III. This object should be also visible 
just above the horizon by the MAGIC telescope in the Ca- 
nary Islands. The expected photon flux above 100 GeV is 
~ 10 -10 cm -2 s -1 which is almost two orders of magnitude 
above the sensitivity of M AGIC in its Large Zenith Angle 
mode ([Blanch et al., 1998|). MAGIC may also probe the ICS 



component of the spectrum (above 0.5 TeV) provided the 
energy of beam particles (i.e. primary electrons) is higher 
than assumed in our model. 

Our ranking list of Figure ^| also includes J0218+4232, 
a pulsar with a likely detection by EGRET at about 3<r level 
(Verbunt et al. 1996, Kuiper & Hermsen 2000). Our theo- 
retical gamma ray flux for this object is a factor of up to ten 
below the expected GLAST sensitivity, mainly because of 
the large distance. Thus within the model presented above 
we remain skeptical to this report about its positive detec- 
tion, knowing that the gamma ray source is coincident with 
an AGN. 

Millisecond pulsars are potentially very interesting 
sources of high energy gamma-rays. The gamma-ray radi- 
ation, if detected, shall provide useful information about 
physical conditions in a pulsar magnetosphere: the primary 
electron energy, the nature of the electric field accelerating 
the primary electrons, and indirectly soft photon density. 
We appeal to include the millisecond pulsars into target lists 
for the upcoming gamma-ray observatories, both space- and 
ground-based. 
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